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Abstract 

We performed two independent counting experiments on a /3-emitting source 
of g2^Sm by measuring tlie 7 photon emitted in a fraction of the decays. For 
counting times ranging from 10~^ to 5.12 x 10^ seconds, our measurements 
show no evidence of deviations from Poissonian behavior and, in particular, 
no sign of 1// noise. These measurements put strong limits on non-Poissonian 
components of the fluctuations for the subset of decays accompanied by 7, 
and corresponding limits for the total number of /? decays. In particular, 
the magnitude of a hypothetical flicker floor is strongly bounded also for 
the (3 decay. This result further constrains theories predicting anomalous 
fluctuations in nuclear decays. 

PACS numbers: 05.40. +j, 02.50.-r, 23.90. +w 
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I. INTRODUCTION 



The statistics of the radioactive decay of heavy nuclei have been the subject of much 
experimental and theoretical work in the past decade. So wide an interest was stimulated 
by the conjecture that, owing to the intrinsic fluctuations of the decay rate, the counting 
statistics could depart from the simple Poissonian behavior 

The experimental results are often conflicting, even for the same kind of source. On the 
one hand, there exist investigations both on a (g|^Am and gi^Po |]12|) and /? decays 



(55 ''Cs [|13[) that confirm the Poissonian nature of these processes. On the other hand, several 
experiments carried out both with a (95^ Am, gl^Pu, and gg^Cm [|I?|-[T^) and with f3 sources 
(g'j'^Tl |T^, 39Y and 3gSr-3gY ||2^), find that the counting variance, for long counting 
periods, is higher than the Poissonian value by more than one order of magnitude. 

This anomalous large variance has been taken as an experimental evidence that the power 
spectrum of the decay-rate fluctuations has a contribution that grows as the inverse of the 
frequency / at low frequencies, in addition to the usual frequency-independent Poissonian 
component. Several mechanisms have been proposed as possible sources of this 1// noise: 
quantum self-interference between the wave packets of the emitted particles solid- 
angle fluctuations and random rearrangements within the source 0, spatial 1// noise in the 
detector |T^. As a matter of fact, the interpretation of the decay experiments reporting a 



variance in excess of the Poisson value is still an open problem 0. 

In previous work pT|P^ we considered the decay statistics of a 7 source (gJ^^Sn). In that 
case, we measured that, for counting periods T longer than one hour, the variance of the 
decay rate significantly deviated from the Poissonian prediction. However, that behaviour 
could be fully explained by taking into account the time dependence of the statistics ^ 
without resorting to any exotic effect |]2l|,|22[|. 



The aim of this article is to extend our experimental study to a different nucleus, 52^ Sm, 
that undergoes (5 decay. There are in fact theoretical claims that deviations from Pois- 
sonian statistics could be caused by self-interference of the emitted particles and that these 
deviations should be present only in (3 decays and not in 7 or o; decays. 



II. EXPERIMENTAL SETUP 

The mean lifetime of g2^Sm is (130 ± 12) years |24|. While most of the nuclei directly 



(3 decay into the ground state of gl^Eu, a small fraction (0.91%) /3-decays into an excited 
state of energy i^exc = (21.532 ± 0.068) keV, which then decays to the ground state (mean 
lifetime: 1.38 x 10~^ seconds). In the 3.45 % of cases, this second fast transition produces 
a 21.532 keV 7-photon: our apparatus has been set up to detect this photon. In summary, 
our measurement has the characteristic of selecting a fraction of the total decays = 
(3. 14 ±0.22) X 10~^): those decays that go through the two-step process, (3 emission followed 
by a 21.532 keV photon [Q. We shall discuss later and in the Appendix why and to what 
extent our results on the statistics of the 7's also carry information on the total statistics of 
the (3 decays. 

The source is a crystal of SmFs containing Sm nuclei (activity: 3.7 GBq) shaped as 
a thin disk (diameter: 14 mm) with an aluminium cap. The aluminium cap, which closes 
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the source, filters out the f3 particles. In the two experiments, which we denote A and B, 
we used the same source at different distances (about 15 and 7 cm, respectively) from the 
detector in order to change the count rate: while Poissonian statistics only depends on the 
total number of counts (rate x time), deviations from the standard case (and/or systematic 
errors) could in principle depend also on the rate (see Ref. |^ for one such an example) 



and it is better to have the possibility of performing this kind of checks. 

In each experiment the photons were detected by a disk-shaped crystal of Nal(Tl) (di- 
ameter: 5 cm) integrally mounted on a photomultiplier tube (PMT): we used a crystal 1 
mm thick in experiment A, and a crystal 2 mm thick in experiment B. 

In both experiments, the output signal from the PMT, after being amplified and shaped 
to a Gaussian pulse, passed through a single- channel analyzer, which selected pulses cor- 
responding to an energy window from 2 to 53 keV. The pulse-shaping time constant was 
0.5/is and the time resolution of the single-channel analyzer 0.6/is. The dead time of the 
entire system was about 2.5/is. The energy window was preliminarily set by means of a 
multi-channel analyzer module; we have verified that no appreciable drift of the window oc- 
curred during the experiments, which lasted 76 days (A) and 19 days (B). We verified that 
the stability of the energy window and of the voltage of the power supply were sufficient to 
keep systematic variations of the counting rate below 0.01%, therefore below the statistical 
fluctuations we measured: only for the longest measurements (total counts of the order of 
10^) the fiuctuations-to-signal ratio was as low as 10~^ (1/vTo^). 

Counting was executed by a programmable multi-channel scaler (MCS) module interfaced 
to an IBM PC, which provides for control and data storage. In experiment A (B), a set of 
40 (38) values of the counting period T was preliminarily defined in the control program 
with T ranging from T^i^ = Ti = 10'^ to T^^^ = T40 = 2^ x 100s (T^^ax = T^g = 2^ x 100s). 
For each value of T the MCS module counted the events occurring in each of 64 consecutive 
periods of length T. Count data were saved on hard disk for further off-line analysis. At 
the end of experiment A (B), data were available as 40 (38) sequences of 64 counts Mj 
{k = 1, . . . , 64), one for each of the prefixed values of T. 



III. RESULTS AND DISCUSSION 



We analyzed the data by computing the average count and the Allan variance (see 
Refs. ^,25] and references therein) as function of the time interval T. All our results originate 



from a single uninterrupted run (for each experiment), have been averaged over the same 
number (64) of consecutive intervals, and are statistically independent (each count has been 
used only once). 

First we verified that the count rate during each experiment had no drifts that could 
bias the Allan variance; in particular, the slow exponential decay of the source could not 
affect the Allan variance at the low count rates we operated |2T|. Therefore, it is consistent 
that we consider a constant average rate. We measured this average rate 

T k=l 
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finding m = (5.3687 ±0.00025) x 10^ count/s in experiment A and m = (2.4262 ±0.00022) x 
lO'^ count/s in experiment B. Being the rate constant, the average count for an interval of 
length T 

1 64 

M{T)^-Y^Ml (2) 

has an expectation value proportional to T: {M{T)) = mT. 

There are two reasons for using the Allan variance, which we estimate with an average 
over 63 consecutive measurements 

1 63 

AiT) ^ ^JIM - (3) 

instead of the usual variance. The first and most important reason is that A{T) is finite 
even when the power spectrum grows as 1// at low frequencies: when non-Poissonian fluc- 
tuations might be present, the Allan variance is then a common choice. We remind that the 
power spectrum of Poissonian fluctuations is independent of frequency: S{f) = 2mT; since 
counts are uncorrelated, the Allan and the usual variance have the same expectation (use 
(MJmJ) oc 6jk and Eq. (|)), namely, the average count: {A(T)) = {M(T)) = mT. How- 
ever, if the fluctuations have a power spectrum S{f) = C// (C is a constant independent 
of /), i.e., we are in presence of 1// noise, the expectation value of A is: {A(T)) = Cln4T^ 
(note the different power of T compared to Poissonian fluctuations), while the usual variance 
is infinite 0. A second additional advantage of using the Allan variance is that it is less 
sensible to drifts of the count rate: the correction is independent of the number of intervals 
(64) and not proportional to it, see the Appendix of Ref. . 

Before discussing our results, we wish to comment on our choice of observing the channel 
of the decay characterized by the emission of a 21.532 keV photon. A more detailed dis- 
cussion can be found in the Appendix. We made this choice because we can control better 
the stability of our measurements when detecting photons than when detecting electrons, 
given our present equipment. However, since one of the motivations of our experiment was 
to study fluctuations in a /3 decay, it is natural to ask to what extent the statistics of the 
7 emission reflects the statistics of the (3 decay. The time delay of the emission is so small 
(mean lifetime of the excited state: 1.38 x 10~^ seconds) compared to the time intervals of 
interest that its effect is negligible. Yet one might worry that the fluctuations of the small 
branching ratio (the fraction of decays that on average emit the 7 is only ^ = 0.000314) 
might overwhelm any exotic effect of the original decay. The explicit calculation reported in 
the Appendix shows that: (1) an upper bound on the flicker floor in the statistics of the 7's 
implies an equal bound on the flicker floor in the statistics of the parent decay; (2) an upper 
bound on the ratio of 1// noise to Poissonian noise in the statistics of the 7's implies a 
corresponding bound for the statistics of the parent decay weaker by a factor 1/,^; (3) upper 
bounds on less singular, e.g. frequency independent, deviations from Poissonian behavior in 
the statistics of the 7's imply corresponding bounds on the parent decay: these bounds on 
the statistics of the P decay are also weaker by a factor 1/^ ~ 3000. 

We report in Figs. |I| and |^ the ratio R{T) = A{T)/M'^{T) (reduced Allan variance) 
versus the inverse of the number of counts 1/M(T) for experiments A and B, respectively. 
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Both experiments show that R{T) depends hnearly on 1/M(T) = l/(mT) with unit slope 
in the range of T considered. 

The data perfectly fit the Poisson prediction R{T) = M/M"^ = 1/M oc l/T; this 
prediction is also reported in Figs. |l| and as a solid line. A fit to the data yields 
M(T) X R{T) = 0.99 ±0.02. 

On the contrary, a power spectrum S{f) = C/ f would yield R(T) = (C In 4 T^)/(mT)^ = 
Cln4/m^. Therefore, if we suppose that both a Poissonian and a 1// contribution are 
present, when T is large enough the Poissonian contribution becomes negligible and R{T) 
goes to a constant [F = Cln4/m^): this constant F is usually called flicker floor. We 
measured values of R{T) as low as 6 x 10^^ (3 x 10^^) in the experiment A (B) at T^ax = 
5.12 X 10*^ s (Tjnax = 1-28 X 10^ s) without seeing deviations from Poissonian behavior and, in 
particular, no signal of the curve turning up at large T and becoming constant. Therefore, 
we conclude that, if a flicker floor is present, F < 3 x 10~^; as discussed in the Appendix, 
this limit is valid also for the f3 decay. 

If we express the power spectrum as the sum of the Poissonian component plus a hy- 
pothetical 1// component: S{f) = 2mT + C/f, in the range of frequencies accessible 
by our experiments (/ > l/Tmax), the limit on the flicker floor implies an upper limit 
on the ratio of the strength of the 1// contribution {C/f) relative to the Poissonian one 
(2mT), i.e., a limit on the ratio (C//)/(2mTmax) < C/(2m): {C / f) / (mT^^^) < 1 x lO'^ 
{C/m < 2.5 X 10^^). These limits on the strength of the 1/f noise are valid for the channel 
of the decay with 7 emission; for the total f3 decays the limit is weaker (see Appendix): 
iCp/f)/{mf3T^,^) < 3 X 10-2. 

The model of quantum 1/f noise proposed by Handel predicts F = 8a( In 2 (Af/c)^/ (37r) 
for (3 decays, see Eq. (3.6) of Ref. |Q, and references therein; here a ~ 1/137 is the fine 
structure constant, < C < 1 is a coherence factor, and Av/c is the velocity change of 
the particles in the emission process relative to the speed of light c: if K/3 is the kinetic 
energy of the electron, (Av/c)^ = 1 — [1 + Kj3/{mc^)]~''^. Since we did not measured the 
electron energy, our data are averaged over the entire electron-energy spectrum. Therefore, 
we can only give an estimate of the limit on the coherence factor by using the average 
electron energy: (K/s) = 13.96 keV. The fact that we do not see any flicker floor implies, 
in the context of Handel's model, that the coherence factor ( must be smaller than about 
10^^. Our limit should be compared to the recent positive determinations of ( in the range 
5.2 X 10"^ < C < 8-3 X 10"^, that Gopala et al. [^] have made, albeit in different f3 decays: 
ggSr, ggY aud 8?^T1. We do not have any explanation why ( should be more than two orders 
of magnitude larger in those decays compared to our upper limit. 

IV. CONCLUSIONS 

We have measured the counting rate of secondary 7 rays from a (3 source of Sm for 
counting periods ranging from 10"^ to 5.12 x 10^ seconds, and studied the fluctuations of 
the rate by means of the Allan variance. 

• We have found no evidence of deviations from Poissonian behavior up to a ratio of 
fluctuations to signal as low as 5 x 10"^. 
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The ratio between a hypothetical 1// component of the power spectrum and the usual 
Poissonian contribution must be less than 1 x 10~^ at the longest time interval (lowest 
frequency) that we have measured (Tmax = 5.12 x 10"^ s). 

We found no evidence of flicker floor. The upper bound on a hypothetical flicker floor 
is 3 X 10~^; this limit is valid also for the statistics of the total (3 decays. 

If our upper bound the flicker floor is interpreted in the context of Handel's theory of 
1// noise predicting coherent interference of the emitted charged particle, the coher- 
ence factor C, for this decay must be less than about 1 x 10~^: this number is more 
than two orders of magnitude smaller than the one that has been recently proposed in 
the literature albeit for different decays pO . 
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APPENDIX: 

In this Appendix, we discuss the relation between fluctuations of the number of total 
decays (in our case, the total number of (3 decays) and fluctuations of the number of decays 
in a subchannel (in our case, the fraction of (3 decays that produce a 7 photon with energy 
21.532 keV). 

The main results of this Appendix are summarized by Eqs. (|A8|) and ([A9|). 

We consider only the effect of the fluctuations of the branching ratio and not the effect of 
time delay between the first and second decay, which in general has the effect of a low-pass 
filter since this second decay (the 7 emission) is practically instantaneous for the case 
under study (mean lifetime 13.8 ns compared to a time resolution of the order of /xs and to 
the shortest time interval considered: 1 ms). 

In the following we shall use the symbol M when referring to the number of detected 
7's and the symbol when referring to the corresponding (total) number of (3 decays. We 
shall also use the subscript 7 (/3) referring to the partial daughter statistics (total parent 
statistics). Let us define two kinds of averages: 

{■ ■ ■)n average over the 7-count distribution keeping the number of j3 counts N fixed; 

( • ■ ■ ) average over the /?-count distribution. 

Using the symbol ^ for the branching ratio, we can write 

{M)m = (Al) 
M={{M)^)=^N , (A2) 

where for simplicity we use the symbol M to indicate the number of 7 counts twice averaged 
both over the 7 and (3 distributions and, at the same time, the symbol to indicate the 
average /3 counts (over the (3 distribution). 
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Since our experiments do not show any deviation from Poissonian behavior, we can 
readily put hmits on non- Poissonian components of the 7 counts. The imphcations for the 
total (3 decay can be assessed by assuming that the 7 distribution at fixed number of 13 
decays N is standard (binomial and frequency independent) and by considering the effect 
of fluctuations of A^: 

((M, - (M)^J(M,. - (M)^J)^^^^ = 5.,(M)^^ = 5,,^{l - OAT, , (A3) 

where the indices i and j indicate different counting intervals. 

We first consider the average at fixed Ni and A'j+i of (Mj — Mj+i)^, which by adding and 
subtracting {M)^. = ^Ni can be written as 

{{M, - M,+i)2);v.7v,+, = ( m - {M)nJ - (M,+i - (M);v,^J + m+i - )n.n.^, (A4) 

= {((M, - {M)^f)^^ + ((M,+i - {M)^^^J')^^^, + e{N, - N,^,y 
-2((M, - (M);vJ(M,+i - (M)^,^J)^,^,^, (A5) 

= e(l-0(iV. + iVm)+e'(iV.-iVm)' , (A6) 

where we have applied Eq. ( [A3D to the first and second line of Eq. ( |A5D , while the third line 
is identically zero. 

If we now divide the above result by two and average it over the (3 distribution, we find 
(considering that for a stationary process (Aj) = A^ independently of i): 

^(((M, - M,+i)')^.^...) = e(l - ON + elUN. - AT.+i)^) . (A7) 



The left-hand side of Eq. ( A7) is the expectation value of the Allan variance of the 7 counts, 
which we measure with the statistics defined in Eq. (^, while the right-hand side is the 
expectation value of the Allan variance of the total /3-decay counts. If we define (Ap) as 
the Allan variance and = / 
(/9) counts, Eq. (|A7|) becomes 



2 2 

the Allan variance and = A.y/M [Rp = Ap/N ) as the relative Allan variance of the 7 



A, = {l-i)M + eAp (A8) 
R,= ^-^ + Rp. (A9) 

which constitute the main result of this Appendix. 

In the following we analyzes the consequences of Eq. (|A9|) for our experimental study. 



a. Flicker floor 

If the (3 decay has a 1// component that produces a flicker floor Fp in the relative Allan 
variance, i.e., Rp = l/N + F/s, the fact that no deviation of R^ from 1/M has been observed 
for R^ as low as 3 x 10~^ implies not only that < 3 x 10~^, but also that Fp < 3 x 10~^, 
where we have used Eq. (^) dropping ^ ^ 3 x 10^^ compared to 1: (1 — ^) ~ 1 and 
1/M + 1/N ^ 1/M. 
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b. 1/f noise 



If wc are interested in the ratio of the 1/f contribution (Cj^/f) relative to the Poissonian 
one (2m^T), we should recall that the rate of the /3 decay is ruf) — m^/^ and the that 
the constant C is related to the flicker floor by C = F x w? /InA. Then this ratio for 
the j3 decay can be related to the same ratio for the 7 decay by using the fact that the 
limit on Fp and the one on are equal: (C/3//)/(2m^Tmax) = {Fprnp/ f) / {T,,,.^^\nA) = 
{1/0 X (F^m-y//)/(Tinaxln4). We loose a factor 1/^ ~ 3000 going from the upper bound on 
the ratio of the 1/f contribution relative to the Poissonian one for the 7 statistics (1 x 10~^) 
to the upper bound on the same ratio for the P statistics (3 x 10~^). 

c. Frequency-independent non-Poissonian component 

If instead we suppose that the f3 decay has a frequency-independent deviation from 
Poissonian statistics, i.e., Rp = k/ N = k^/M, the fact that no deviation of i?^ = (1 — ^ + 
O/M from 1/M has been observed (M x = 0.99 ± 0.02) implies also that (1 - k)^ = 
0.01 ± 0.02 and, consequently, that k, — —30 ± 60. 

In conclusion, we have show in this Appendix that a measurement of a process (/3 decay) 
by selecting a subprocess (detecting the 7 emitted in a fraction of the decays) whose branch- 
ing ratio ^ is itself a statistical variable corresponds, as might have been expected, to the 
use of a detector with efficiency not greater than ^. Therefore, we loose a factor 1/^ in most 
limits on dimensionless quantities when passing from statistics of the subprocess to total 
statistics of entire process. However, there exist quantities, such as the flicker floor, that can 
be determined from the partial statistics without loosing any sensibility. The reason of this 
different behavior is related on how strongly the noise under study depends on the number 
of events N compared to the usual y/N dependence. 
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FIGURES 

FIG. 1. Relative Allan variance R{T) versus the inverse mean number of decays 1/M (lower 
scale) and the inverse time interval 1/T (upper scale) for the experiment A (average rate 
m = M/T = 5.3687 x 10^). Diamonds are the experimental values. The solid line is the Poissonian 
prediction R = 1/M. 

FIG. 2. Same as Fig. Q for the experiment B (average rate m = 2.4262 • 10^). 
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